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Abstract 

We validate the baryonic Tully Fisher (BTF) relation by exploring the Tully Fisher (TF) and BTF properties 
of optically and Hl-selected disk galaxies. The data includes galaxies from: Sakai et al. (2000) calibrator sample; 
McGaugh et al. (2000: MC2000) I-band sample; and 18 newly acquired Hl-selected field dwarf galaxies observed 
with the ANU 2.3 m telescope and the ATNF Parkes telescope from Gurovich's thesis sample (2005). 

As in MC2000, we re-cast the TF and BTF relations as relationships between baryon mass and W2q- First 
we report some numerical errors in MC2000. Then, we calculate weighted bi-variate linear fits to the data, and 
finally we compare the fits of the intrinsically fainter dwarfs with the brighter galaxies of Sakai et al. (2000). 
With regards to the local calibrator disk galaxies of Sakai et al. (2000), our results suggest that the BTF 
relation is indeed tighter than the TF relation and that the slopes of the BTF relations are statistically flatter 
than the equivalent TF relations. Further, for the fainter galaxies which include the I-band MCG2000 and 
Hl-selected galaxies of Gurovich's thesis sample, we calculate a break from a simple power law model because 
of what appears to be real cosmic scatter. Not withstanding this point, the BTF models are marginally better 
models than the equivalent TF ones with slightly smaller x?ed- 
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1 Introduction 

McGaugh et al. (2000: MC2000) provides a valuable sam- 
ple of I-band data for 63 optically selected dwarf galaxies. 
Several of these galaxies are isolated field dwarfs and a 
few are among the slowest rotators of the known popula- 
tion of disk galaxies. Using these galaxies, supplemented 
by other brighter galaxies, MC2000 conjected that a log 
vs. log baryonic TF relation (BTF) which includes all the 
visible baryons (gas plus stars) appears to be more nearly 
linear than the TF relation which includes only the stellar 
component. If true, a BTF relation relates the baryon con- 
tent to the dark matter halo's of galaxies, and was probably 
established at the time of their formation. 

Gurovich, Freeman, and McGaugh (in preparation) 
have shown that arithmetic errors crept into the MC2000 
I-band data. The I-band magnitudes were systematically 
over- luminous by 2 x (1/ — /) (sign error) and the HI masses 
were overestimated by 0.11 dex {Hq error). In order to 
establish the validity of the BTF relation we went back 
to the original published data and reanalysed it following 
the methodology of MC2000. The empirical BTF relation 
(Freeman, 1999; MC2000) for disk galaxies relates the to- 
tal baryon mass (gas + stars) and VK20, the 20th percentile 
width of the HI profile. For the brighter galaxies W20 is a 
fairly direct measure of the amplitude of the flat part of 
the rotation curve, but for the less massive and more slowly 
rotating galaxies this is not so clear. 
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2 Calculating Stellar and Gas Baryons 

We went to the source of the MCG2000 I-band data for 
the 64 dwarf galaxies: Pildis et al. (1997), Schombert et 
al. (1997) and Eder & Schombert (2000), and, from appar- 
ent magnitudes kindly made available by James Schombert 
(private communication), we recalculated the stellar and 
gas baryon masses following the MCG2000 formulation. 
For all 64 galaxies we measured the velocity widths W20 
based on the spectra of Eder & Schombert (2000) and cor- 
rected for the incUnations by Pildis et al. (1997). 

The MC2000 data come from an optically selected 
sample of dwarf galaxies. In this analysis we have also 
included 18 Hl-selected field dwarf galaxies with inclina- 
tions > 45° from the HIPASS survey, which are part of 
a larger thesis study by Gurovich (2005). In this paper 
we have also included the brighter local calibrator sam- 
ple (with Cepheid distances) of Sakai et al. (2000). For 
the Hl-selected galaxies we collected H-band photometry 
from the Siding Spring 2.3 m telescope and obtained VK20 
and HI masses from narrow-band (8 MHz, 1024 channels) 
observations with the Parkes telescope. All distances were 
determined using cosmological redshifts, from the HI sys- 
temic velocities, calculated for Hq — 75 kms~^Mpc~^, ex- 
cept where reliable secondary distances were available. The 
HI data for the Sakai et al. (2000) galaxies were calcu- 
lated from Martin (1998) and Pilyugin et al. (2004) but 
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Figure 1. TF (left panel) and BTF (right panel) plots for the, corrected I-band MC2000 and adapted Sakai et al. (2000) I-band 
calibrator, data showing the correlation between stellar mass and velocity width W20- Stellar masses are represented by circles. "Total" 
masses, which include the stellar plus gas components, are represented by diamonds. The stellar and "total" mass bi-variate weighted 
fits are shown. It is clear that the BTF fits are fiatter but, that the TF and BTF relations break for the fainter MCG2000 sample. 
For the TF and BTF Sakai I-band data the x^ed ^-^^^ ^-^^^ ^^^^ slopes of 3.7 ± 0.3 and 3.0 ± 0.3 respectively and for the 
TF and BTF MCG2000 I-band data the x^ed '^^^ ^^^^ slopes of 6.4 ± 0.3 and 5.5 ± 0.3 respectively. 



# Sakai et al. (2000) H-band stellar baryons 

— Sakai et al. (2000) weighted fit 

# Gurovich PhD thesis sample H-band stellar baryons 

— Gurovich PhD thesis sample weighted fit 



^ modified Sakai et al. (2000) H-band stellar + gas baryons 
— modified Sakai et al. (2000) weigthed fit 
<^ Gurovich PhD thesis sample H-band stellar -i- gas baryons 
Gurovich PhD thesis sample weighted fit 



= (M/L), 



0:H ^©:H 



: 0.6 * L 



'©:H 



"©:gas 



■(M/L)^^„*L^^„= 



1.4*M^.j^^-F 0.6 ' 



11 - 



10 - 



© 



O 



9 - 



I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



7 - 




I I I I I I 




O 



- 7 



111/11 I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 1 1 h I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 

1.6 1.8 2 2.2 2.4 2.6 2.8 3 1.6 1.8 2 2.2 2.4 2.6 2.8 3 

log (km/s) log (km/s) 

Figure 2. TF (left panel) and BTF (right panel) plots for the H-band data of the Hl-selected galaxies from Gurovich's thesis sample 
and, for the adapted Sakai et al. (2000) H-band calibrator data. Stellar masses are represented by circles. "Total" masses, which 
include the stellar plus gas components, are represented by diamonds. The stellar and "total" mass bi-variate weighted fits are shown. 
It is clear that the slopes of the BTF fits are fiatter, but that the TF and BTF relations breakdown for the fainter Hl-selected sample. 
For the Hl-selected sample the TF and BTF x^ed ^•'^ with slopes 5.5 ± 0.9 and 2.2 ± 0.4 respectively and for the Sakai 

H-band data the TF and BTF Xred ^-^^^ ^-^^^ "^^^^ slopes of 4.4 ± 0.3 and 3.3 ± 0.3 respectively. 
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Figure 3. Gas fraction for the optically selected samples of 
MC2000 and Sakai et al. (2000) and Hl-selected sample of galax- 
ies in this paper. 

rescaled to the Sakai et al. (2000) distances. Total magni- 
tudes are extrapolated to infinity. We performed a TF and 
BTF analysis on all four data sets by following MC2000 
and adopted M / L\ — 1.7, M / L^i — 0.6 and estimated the 
gas mass Mg = 1.4Mhi (from cosmological abundances). 

2.1 Weighted Bi-variate Fits 

We used f itexy in Numerical Recipes (Press et al. 1992) to 
find weighted bi-variate least-squares linear TF and BTF 
fits. This algorithm fits the line that minimizes the scatter 
in both the log(VK2o) and the logM variables, weighting 
the residuals by the uncertainties in each quantity. 

2.2 Errors Used for Weighting the Bi-variate 
Fits 

All estimates of baryon masses and have errors arising 
from uncertainties in the: assumed distances; VK20 measure- 
ments; integrated HI fluxes (momO); inclinations and total 
extrapolated magnitude measurements. 

For the Gurovich sample and the corrected MCG2000 
I-band data, we adopt random peculiar velocities of 
180 kms~^ (2.4 Mpc in distance) from Tonry et 
al. (2000), integrated momO errors of 20% and to- 
tal extrapolated magnitude errors of ~ 0.1 mag. We 
believe that our distance errors are probably over- 
estimates (Karachentsev et al. 2003). The er- 
rors are in the range Skms""*^ <^ 5W2q ^ 25 km s""*^ and 
8kms~"^ < SW20 ^ 50kms~"^ for the two samples, respec- 
tively. The errors for the Sakai et al. (2000) galaxies were 
taken from that paper and the integrated HI flux errors 
from the standard deviation of the mean values found in 
the literature re-scaled to the Sakai et al. (2000) distances. 
These errors were propagated following Bevington (1969) 
and subsequently used in the weighted bi-variate fits. 

3 Results of the Fits and Conclusions 

The fits are plotted in Figs. Q] and |21 The slope and x?ed 
can be found in the captions. These results indicate that for 
the brighter spiral galaxies (Sakai et al. 2000), the BTF fits 
are tighter than the TF fits and have shallower slopes than 
the equivalent TF fits. The TF and BTF fits for the fainter 
galaxies of MCG2000 and of the Gurovich thesis sample 
break from this simple power law model as is evident from 
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the large xLd; X?ed,TF = 4.4 and x?ed,BTF =4.1 for the I- 
band and x?ed,TF — 9-7 and X?ed,BTF — 6.9 for the H-band 
data, respectively. We would like to highlight the fact that 
careful consideration went into the error analysis and that 
this result reveals an inadequate model and is not due to 
an underestimation of the errors. 

4 Discussion 

Figure 3 shows the gas mass fraction (G/S) for the galaxies 
in this paper. As is expected, the fainter galaxies are rel- 
atively more gas-rich than the brighter spirals of Sakai et 
al. (2000). It is interesting that the brighter spirals which 
have smaller (G/S) are better modeled by a BTF than 
a TF model. This is consistent with the conjecture that 
there is a strict proportionality between the HI and dark 
Baryons: MCG2000 and Pfenniger and Revaz (2004). It 
is equally interesting that the fainter disk galaxies that 
have larger (G/S) are not well modelled by a simple BTF 
power law fit. It could be that this non-linearity arises 
partly from the assumptions of the constant M/Li — 1.7, 
M/Lh — 0.6 and Mg - IAMhl In fact Pfenniger and Re- 
vaz (2004) conclude that a Mg = 3Mhi is a more realistic 
proportionality. It is perhaps worth discussing the fainter 
galaxies in the right panels of Figs. 1 and 2 that appear to 
be exceptionally baryon-rich for their observed W2o- 

Are these galaxies unusual in having a large baryon 
mass or a small value of VK20? One possible interpretation 
is that the gas in these galaxies is not quite at the flat 
part of the rotation curve and hence their W20 values are 
lower limits of the true rotation amplitude. However, we 
note that these galaxies do appear to better follow the 
(stellar) TF relation. Another possibility is that the gas 
fell into these galaxies after the stellar disk had settled, so 
that the extra baryons are not part of the original galaxy. 
An even more speculative possibility is that a subsample 
of these systems have dark halos which are significantly 
less dense than would be expected from the Kormendy & 
Freeman (2004) scaling laws for dark halos, and hence have 
atypically low values of W20 ■ 
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